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Abstract In this paper, we study the distribution and behaviour of internal 
equilibria in a d-player n-strategy random evolutionary game where the game 
payoff matrix is generated from normal distributions. The study of this paper 
reveals and exploits interesting connections between evolutionary game theory 
and random polynomial theory. The main contributions of the paper are some 
qualitative and quantitative results on the expected density, fn,d, and the 
expected number, E{n, d), of (stable) internal equilibria. Firstly, we show that 
in multi-player two-strategy games, they behave asymptotically as -\/c? — 1 as d 
is sufficiently large. Secondly, we prove that they are monotone functions of d. 
We also make a conjecture for games with more than two strategies. Thirdly, 
we provide numerical simulations for our analytical results and to support the 
conjecture. As consequences of our analysis, some qualitative and quantitative 
results on the distribution of zeros of a random Bernstein polynomial are also 
obtained. 

Keywords Random evolutionary games • Internal equilibria • Random 
polynomials • Multi-player games 


1 Introduction 

l.I Motivation 

Evolutionary game theory (EGT) has been proven to be a suitable mathemat¬ 
ical framework to model biological and social evolution whenever the success 
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of an individual depends on the presence or absence of other strategies [34, 
27,39]. EGT was introduced in 1973 by Smith and Price [35] as an applica¬ 
tion of classical game theory to biological contexts, and has since then been 
widely and successfully applied to various fields, not only biology itself, but 
also ecology, population genetics, and computational and social sciences [34, 
4,27,39,10,41,23,17]. In these contexts, the payoff obtained from game inter¬ 
actions is translated into reproductive fitness or social success [27,39]. Those 
strategies that achieve higher fitness or are more successful, on average, are 
favored by natural selection, thereby increase in their frequency. Equilibrium 
points of such a dynamical system are the compositions of strategy frequencies 
where all the strategies have the same average fitness. Biologically, they pre¬ 
dict the co-existence of different types in a population and the maintenance 
of polymorphism. 

As in classical game theory with the dominant concept of Nash equilibrium 
[37,36], the analysis of equilibrium points in random evolutionary games is of 
great importance because it allows one to describe various generic properties, 
such as the overall complexity of interactions and the average behaviours, 
in a dynamical system. Understanding properties of equilibrium points in a 
concrete system is important, but what if the system itself is not fixed or 
undefined? Analysis of random games is insightful for such scenarios. To this 
end, it is ambitious and desirable to answer the following general questions: 

How are the equilibrium points distributed? How do they behave 
when the number of players and strategies change? 

Mathematical analysis of equilibrium points and their stability in a gen¬ 
eral (multi-player multi-strategy) evolutionary game is challenging because one 
would need to cope with systems of multivariate polynomial equations of high 
degrees (see Section 2 for more details). Nevertheless, some recent attempts, 
both through numerical and analytical approaches, have been made. One ap¬ 
proach is to study the probabilities of having a concrete number of equilibria, 
whether all equilibrium points or only the stable ones are counted, if the pay¬ 
off entries follow a certain probability distribution [16,22]. This approach has 
the advantage that these probabilities provide elaborate information on the 
distribution of the equilibria. However, it consists of sampling and solving of 
a system of multivariate polynomial equations; hence is restricted, even when 
using numerical simulations, to games of a small number of players and/or 
small number of strategies: it is known that it is impossible to (analytically) 
solve an algebraic equation of a degree greater than 5 [1]. Another possibil¬ 
ity is to analyze the attainability of the patterns and the maximal number 
of evolutionarily stable strategies (ESS) [34,27], revealing to some extent the 
complexity of the interactions. This line of research has been paid much atten¬ 
tion in evolutionary game theory and other biological fields such as population 
genetics [28,48,20,7,9,2,16,22,17]. More recently, in [14], the authors investi¬ 
gate the expected number of internal equilibria in a multi-player multi-strategy 
random evolutionary game where the game payoff matrix is generated from 
normal distributions. By connecting EGT and random polynomial theory, they 
describe a computationally implementable formula of the mean number of in- 
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ternal equilibrium points for the general case, lower and upper bounds the 
multi-player two-strategy random games, and a close-form formula for the 
two-player multi-strategy games. 

In this paper, we address the aforementioned questions, i.e., of analysing 
distributions and behaviours of the internal equilibria of a random evolutionary 
game, in an average manner. More specifically, we first analyse the expected 
density of internal equilibrium points, fn.di't), i.e. the expected number of such 
equilibrium points per unit length at point t, in a d-player n-strategy random 
evolutionary game where the game payoff matrix is generated from a normal 
distribution (for short, normal evolutionary games). Here the parameter t = 
denotes the ratio of frequency of strategy i S {1, • • • , n — 
1} to that of strategy n, respectively (more details in Section 2). In such a 
random game, we then analyse the expected number of internal equilibria, 
E(ji,d), and, as a result, characterize the expected number of internal stable 
equilibria, SE{n,d). We obtain both quantitative (asymptotic formula) and 
qualitative (monotone properties) results of fn,d and E{n,d), as functions of 
the ratios, t, the number of players, d, and that of strategies, n. 

To obtain these results, we develop further the connection between EGT 
and random polynomial theory explored in [14] , and more importantly, estab¬ 
lish appealing (previously unexplored) connections to the well-known classes 
of polynomials, the Bernstein polynomials and Legendre polynomials. In con¬ 
trast to the direct approach used in [16,22], our approach avoids sampling and 
solving a system of multivariate polynomial equations, thereby enabling us to 
study games with large numbers of players and/or strategies. 

We now summarise the main results of the present paper. 


1.2 Main results 

The main analytical results of the present paper can be summarized into three 
categories: asymptotic behaviour of the density function and the expected 
number of (stable) equilibria, a connection between the density function with 
the Legendre polynomials, and monotonic behaviour of the density function. 
In addition, we provide numerical results and illustration for the general games 
when both the numbers of players and strategies are large. 

To precisely describe our main results, we introduce the following notation 
regarding asymptotic behaviour of two given functions u and v : [0, -l-oo) —>■ 
[0, -l-oo) 

^ T there exists a positive constant C such that for all k S [0, -l-oo) 

u{k) < Cv{k), 

u ~ u it holds that u<v and v <u. 

Note that throughout the paper we sometimes put arguments of a function 
as subscripts. For instance, the expected density of internal equilibrium points. 
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/n,d(t), besides t, is also analyzed as a function of n and d. We will explic¬ 
itly state which parameter(s) is being varied whenever necessary to avoid the 
confusion. 

The main results of the present paper are the following. As described above, 
fn,d{'t) denotes the expected number of internal equilibrium points per unit 
length at point t, in a d-player n-strategy random evolutionary game where 
the game payoff matrix is generated from a normal distribution; E(n,d) the 
expected number of internal equilibria; and SE{n, d) the expected number of 
internal stable equilibria. The formal dehnitions of these three functions are 
given in Section 2. 

(1) In Theorem 1, we prove the following asymptotic behaviour of f 2 ,d{t) for 

all t > 0: f 2 ,d{t) ~ y/d — 1- We also prove that f 2 ,d{i) is always bounded 
from above and lim^-^oo = 0- 

(2) In Theorem 2, we prove a novel upper bound for the expected number 

of multi-player two-strategy random games, E{2,d) < ^/d — 1 ln(d — 1) 
and obtain its limiting behaviour: lim This upper bound 

is sharper than the one obtained in [14, Theorem 2], which is, E{2,d) < 
(d — 1)3. These results lead to two important corollaries. First, we obtain 
a sharper bound for the expected number of stable equilibria, SE(2, d) < 
I\/(i — 1 ln(d — 1), and the corresponding limit, lim f see 

^ d—>-oo ^ 

Corollary 1. The second corollary. Corollary 2, is mathematically signifi¬ 
cant, in which we obtain lower and upper bounds and a limiting behaviour 
of the expected number of real zeros of a random Bernstein polynomial. 

(3) In Theorem 3, we establish an expression of f 2 ,d{t) in terms of the Legendre 
polynomial and its derivative. 

(4) In Theorem 4, we express f 2 ,d{t) in terms of the Legendre polynomials of 
two consecutive order. 

(5) In Theorem 5, we prove that is a decreasing function of d for any 

given t > 0. Consequently, and are decreasing functions of 

d. 

(6) In Proposition 2, we provide a condition for f 2 ,d{t) being an increasing 
function of d for any given t > 0. We conjecture that this condition holds 
true and support it by numerical simulation. 

(7) In Theorem 6, we provide an upper bound for /n, 2 (t)- We also make a 
conjecture for /n,d(t) and E{n,d) in the general case {n,d > 3). 

(8) We offer numerical illustration for our main results in Section 4.2. 

The density function fn,d(d‘) provides insightful information on the distribution 
of the internal equilibria: integrating fn.di't) over any interval produces the 
expected number of real equilibria on that interval. In particular, the expected 
number of internal equilibria E{n, d) is obtained by integrating /„,d(t) over the 
positive half of the space. Theorem 5 and Proposition 2, which are deduced 
from Theorems 3 and 4, are qualitative statements, which tell us how the 
expected number of internal equilibria per unit length f 2 ^d in a d-player two- 
strategy game changes when the number of players d increases. Theorem 1 
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quantifies its behaviour showing that f 2 ^d is approximately (up to a constant 
factor) equal to ^/d — 1. The function h{d) := \/d — 1, as seen in Theorem 1, 
certainly satisfies the properties that h{d) increases but decreases. Thus, 
it strengthens Theorem 5 and further supports Conjecture 1. Theorem 2 is 
also a quantitative statement which provides an asymptotic estimate for the 
expected number of internal (stable) equilibria. 

Furthermore, it is important to note that the expected number of real zeros 
of a random polynomial has been extensively studied, dating back to 1932 with 
Block and Polya’s seminal paper [6] (see, for instance, [15] for a nice exposition 
and [45,38] for the most recent progress). Therefore, our results, in Theorems 
2, 3 and 4, provide important, novel insights within the theory of random 
polynomials, but also reveal its intriguing connections and applications to 
EGT. 


1.3 Organisation of the paper 

The rest of the paper is structured as follows. In Section 2, we recall relevant 
details on EGT and random polynomial theory. Section 3 presents full analysis 
of the expected density function and the expect number of internal equilibria 
of a multi-player two-strategy game. The results on asymptotic behaviour and 
on the connection to Legendre polynomials and its applications are given in 
Sections 3.1 and 3.2, respectively. In Section 4, we provide analytical results for 
the two-player multi-strategy game and numerical simulations for the general 
case. Therein we also make a conjecture about an asymptotic formula for the 
density and the expected number of internal equilibria in the general case. In 
Section 5, we sum up and provide future perspectives. Finally, some detailed 
proofs are presented in the appendix. 


2 Preliminaries: replicator dynamics and random polynomials 

This section describes some relevant details of the EGT and random polyno¬ 
mial theory, to the extent we need here. Both are classical but the idea of 
using the latter to study the former has only been pointed out in our recent 
paper [14]. 


2.1 Replicator dynamics 

The classical approach to evolutionary games is replicator dynamics [46,50, 
27,44,39], describing that whenever a strategy has a fitness larger than the 
average fitness of the population, it is expected to spread. Formally, let us 
consider an infinitely large population with n strategies, numerated from 1 to 
n. They have frequencies a;^, 1 < i < n, respectively, satisfying that 0 < 3:^ < 1 
and X]r=i interaction of the individuals in the population is in 
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randomly selected groups of d participants, that is, they play and obtain their 
fitness from d-player games. We consider here symmetrical games (e.g. the 
public goods games and their generalizations [24,26,43,40,21]) in which the 
order of the participants is irrelevant. Let be the payoff of the focal 

player, where io (1 < io < n) is the strategy of the focal player, and ik (with 
1 < and 1 < fc < d — 1) be the strategy of the player in position k. 

These payoffs form a (d — l)-dimensional payoff matrix [16], which satisfies 
(because of the game symmetry) 




( 1 ) 


whenever {i[ ..., is a permutation of {zi..., id-i}- This means that only 
the fraction of each strategy in the game matters. 

The equilibrium points of the system are given by the points (xi,..., a;„) 
satisfying the condition that the fitnesses of all strategies are the same, which 
can be simplified to the following system of n — 1 polynomials of degree d — 1 
[14] 


E 



ki—d—1 

i = l 



for i = 1,..., n — 1, (2) 


where := - otl and (\ 

k = l 

multinomial coefficients. Assuming that all the payoff entries have the same 
probability distribution, then all * = 1, ■ • ■, — 1, have symmetric 

distributions, i.e. with mean 0 [22]. 

We focus on internal equilibrium points [16,22,14], i.e. 0 < Xi < 1 for all 
\ < i < n — 1. Hence, by using the transformation ti = with 0 < ti < +oo 
and 1 < f < n — 1, dividing the left hand side of the above equation by x^“^ 
we obtain the following equation in terms of (ti,... ,tn-i) that is equivalent 
to (2) 



TJ.-1 

ki<d—l 

( 3 ) 

Hence, finding an internal equilibria in a general n-strategy d-player random 
evolutionary game is equivalent to find a solution (yi,... ,y„_i) S of 

the system of {n—1) polynomials of degree (d—1) in (3). This observation links 
the study of generic properties of equilibrium points in EGT to the theory of 
random polynomials. 
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2.2 Random polynomial theory 


Suppose that all are Gaussian distributions with mean 0 and vari¬ 

ance 1, then for each i {\ < i < n — 1), 


multivariate normal random vector with mean zero and covariance matrix C 
given by 

2 \ 


C = diag 



(4) 


0<fci<d-l, Y. ki<d-l 


The density function fn^d and the expected number E(ri, d) of equilibria can 
be computed explicitly. The lemma below is a direct consequence of [15, 
Theorem 7.1] (see also [14, Lemma 1]). For a clarity of notation, we use 
bold font to denote an element in high-dimensional Euclidean space such as 

t = (ti, . . . ,tn_l) e 


Lemma 1 Assume that are independent normal random vectors 

with mean zero and covariance matrix C as in (4). The expected density of 
real zeros of Eq. (3) at a point t = (ti,..., tn-i) is given by 


fn,d(t)='^ "4" (^) (detL(t)) 2 , (5) 


where E denotes the Gamma function and L(t) the matrix with entries 


Lij(t) 


^(log.(x)rc.(y)) 


y=x=t’ 


with 


viKfCviy) 



n-1 

ki<d—l 

i = l 


( 6 ) 


As a consequence, the expected number of internal equilibria in a d-player 
n-strategy random game is determined by 


E{n,d)= [ fn,d{t)dt. (7) 

Note that the assumption in Lemma 1 is quite limited when applying to games 
with more than two strategies as in that case the independence of the a terms 
does not carry over into the independence of /? terms, see Remark 3 for a 
detailed discussion. 
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3 Multi-player two-strategy games 

We provide mathematical analysis of the expected density function f 2 ,d{t) and 
the expected number of equilibria E{2, d) for a multi-player two-strategy game. 
Section 3.1 presents asymptotic behaviour. A connection to Legendre polyno¬ 
mials and its applications are given in 3.2. In Section 3.2, further applications 
of this connection to study monotonicity of the density function are explored. 


3.1 Asymptotic behaviour of the density and the expected number of 
equilibria 


In the case of multi-player two-strategy games (n = 2), the system of polyno¬ 
mial equations in (2) becomes a univariate polynomial equation 


d-l 




y"(l-2/)“ = 0, 


( 8 ) 


where y is the fraction of strategy 1 (i.e., 1 — y is that of strategy 2) and 
(3k is the payoff to strategy 1 minus that to strategy 2 obtained in a d-player 
interaction with k other participants using strategy 1. It is worth noticing that 

Vd-fc = o,...,d-i, (9) 


is the Bernstein basis polynomials of degree d—1 [42,32]. Therefore, the left- 
hand side of (8) is a random Bernstein polynomial of degree d—1. As a 
by-product of our analysis, see Theorem 2, we will later obtain an asymptotic 
formula of the expected real zeros of a random Bernstein polynomial. 

Letting t = {t G K+), Eq. (8) is simplified to 


d-l 





( 10 ) 


The expected density of real zeros of this equation at a point t is f 2 ,d(t)- 
For simplicity of notation, from now on we write /d(t) instead of f 2 ,d(t). We 
recall some properties of the density function /d(t) from [14, Proposition 1] 
that will be used in the sequel. 


Proposition! ([14]) The following properties hold 

1) fd{0) = ^, = ^ 72 ^ 3 - 

2) fd{t) = ^ where 


G{t) = t 


Md{t) 



log Md{t) 


with 


d-l 

Md(t) = ^ 
fe =0 




( 11 ) 
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3) fd{t) has an alternative representation 

(2-/.(t)f = (12) 

i—1 ^ 
d-1 

where ri > 0 satisfies that Md{f) = 11 

2—1 

f) t fd{t) is a decreasing function. 

5) fd{\) =t^f{t). 

6) E{2, d) = /~ fait) dt = 2 /q' fd{t) dt. 

Example 1 Below are some examples of fd(t), with c? = 2,3 and 4 

^ ,,, 2 VT+f^T¥ ^ 3 Vl + 4t2 + 10^4 + 4^6 + t8 

“ 7rl + t2’ ^ 1 + 4^2 i + gt2+gt^ + t^ 

We recall that t = where y is the fraction of strategy 1. We can write 

the density in terms of y using the change of variable formula as follows. 

Define gd{y) to be the density on the right-hand side of the above expression, 
i.e., 

( 13 ) 

The following lemma is an interesting property of the density function 
which explains the symmetry of the strategies (swapping the index labels con¬ 
verts an equilibrium at y to one at \ — y). Numerical simulations in Section 
4.2 further illustrate this property (see Figure 6). 

Lemma 2 The function y i—> gd{y) is symmetric about the line y = ^, i.e., 

gd{y) = 5d(l - y)- 


Proof We have 

s4i-!/) = A(t((^)4=/,(^) 


1 


- (1 - 


- y) {1 - yy 


where we have used the fifth property in Proposition 1 to obtain the second 
equality above. 


The following theorem provides an upper bound and asymptotic behaviour 
for fd{t) and its scaling with respect to d. 
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Theorem 1 (Asymptotic behaviour of the density function) The fol¬ 
lowing statement holds 


fd{t) < min 


^/d— 1 d—1 


for t > 0. 


(14) 


27rt TT 
As a consequence, for any given t > 0 

h 

d- 

and furthermore 


Proof It follows from the first and the fourth properties in Proposition 1 that 

d-1 


1 

< , 

TT 

y fd{t) „ 

lim , =0, 

d^oo Cl — 1 

(15) 

fd(t) ' 

Vd — 1 . 

(16) 


fdit) < fM = 


(17) 


On the other hand, according to the third property in Proposition 1, we have 


Since it follows that for any t 0 

, , ,,9 1 ^1 

(27r/d(t)) <^- = —, 

2=1 


which is equivalent to fd{t) < ■ Hence, the upper bound of fd{t) in (14) 

holds. 

The upper bound and limit in (15) are obvious consequences of (14). It 
remains to prove the asymptotic behaviour of fd{t) as a function of d. Using 
(14) and the fourth property in Proposition 1, it follows that, for 0 < t < 1 


d-1 1 

27r 


/d(l) < fd{t) < 


y/dTAl 

2711 ’ 


and for t > 1 


1 d- 1 1 

f^ 27r v^2d - 3 


^ ^ 2 /^ (i) 


fd{t) < /d(l) 


d-1 1 

27r V2d^' 


From these estimates, we deduce that fd{t) ^ y'd— 1 for any t > 0. 
We numerically illustrate Theorem 1 in Section 4.2, see Figures 1. 
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Theorem 2 (Asymptotic behaviour of E{2,d)) It holds that 

Vd^<E{2,d)<Vd^\n{d-l). (18) 


Furthermore, we obtain the following asymptotic formula for E(2, d) 

\nE{2,d) _ 1 
rf—fcjo ln(fi — 1) 2 


(19) 


We first provide a proof of this theorem. Then we discuss its implications 
for the expected number of stable equilibrium points in a random game and 
that of real zeros of a random Bernstein polynomial in Corollaries 1 and 2. A 
comparision with well-known results in the theory of random polynomials is 
given in Remark 1. 

Proof The lower bound was derived previously in [14, Theorem 2]. For the 
sake of completeness, we provide it again here. Using the sixth and the fourth 
properties in Proposition 1, we have E{2,d) = 2 fd{t) dt > 2 f{l)dt = 
2/(1) = Therefore, ^d^ < E{2,d). 

The underlying idea of the proof for the upper bound is to split the integral 
range in the formula of E(2,d) into two intervals. The first one is from 0 to 
a, for some a S (0,1); we then estimate fd{t) in this interval by /d(0). The 
second one is from a to 1, which is estimated using the upper bound of fd{t) 
given in (14). The value of a will then be optimized. 


r nd 

E{2,d) = 2 fd{t)dt = 2 / fd{t)dt+ fd{t)dt 

*/ 0 J 0 J Oi 


/O 

< 2 
1 r 


a«0) + ^ / \dt 


2{d — 1) a — Vd — 1 In < 


( 20 ) 


Let h{a) be the expression inside the square brackets in the right-hand side 
of (20). To obtain the optimal (i.e. smallest) upper bound, we minimize h{a) 
with respect to a. The optimal value of a satisfies the following equation 


^h{a) = 2{d - 1) - = 0, 

da a 


which leads to a = 



Substituting this value into (20), we obtain 


E{2,d) < 


^/d^ 


1 -b ln2 -I- - ln(d — 1) 


It follows that E{2,d) < y/d — 1 ln(d — 1), which is (18). 

We now prove (19). By taking logarithm in (18), we obtain 

- ln(d — 1) < hvE{2, d) < - ln(d — 1) -I- lnln(d — 1). 
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Therefore 

1 ^ ln£’(2,(i) ^ 1 lnln(c?—1) 

2 ^ ln(d- 1) 2 ln(d - 1) 

Since Jim = 0, we achieve (19). 


( 21 ) 


Theorem 2 has two interesting implications about the expected number of 
stable equilibrium points in a random game and that of real zeros of a random 
Bernstein polynomial. 


Corollary 1 The expected number of stable equilibrium points in a random 
game with d players and two strategies, SE{2,d), satisfies 


^ SE{2,d) < ^Vd- lln(d- 1), 
and furthermore, satisfies the following limiting behaviour 


\iiSE(2,d) 1 

lim ———-— = 

d^oo \n[d — 1 ) 2 


( 22 ) 


(23) 


Proof From [22, Theorem 3], it is known that an equilibrium in a random game 
with two strategies and an arbitrary number of players, is stable with prob¬ 
ability 1/2. Thus, SE{2,d) = Hence, the corollary is clearly followed 

from Theorem 2. 


Corollary 2 The expected number of real zeros, Eq, of a random Bernstein 
polynomial 

fc ~ 

k^O ^ ^ 

where fit independent standard normal distributions, satisfies 

Vd-l<EB<Vd - lln(d- 1), lim = 7 ^- (24) 

d—>^oo ln(^a — 1) 2 

Proof As mentioned beneath (8), solving B{x) = 0 is equivalent to solving 
Eq. (10). It follows that fd{t) is the expected density of real zeros of B{x). 
Therefore, Eq given by 


/ OO POO 

fd{t)dt = 2 fd{t)dt = 2E{2,d). (25) 

-OO d 0 

Note that the second equality in (25) holds because fd(t) is even in t due 
to (12). The asymptotic behaviour (24) of Ag is then followed directly from 
Theorem 2. 














Title Suppressed Due to Excessive Length 


13 


Remark 1 The study of the distribution and expected number of real zeros of 
a random polynomial is an active research field with a long history dating back 
to 1932 with Block and Polya [6], see for instance [15] for a nice exposition 
and [45,38] for the most recent results and discussions. Consider a random 
polynomial of the type 

d-l 


Vdiz) = ^Ci^iZ\ 


(26) 


i=0 


The most well-known classes of polynomials studied extensively in the litera¬ 
ture are: flat polynomials or Weyl polynomials for elliptic polynomials 


(EP) or binomial polynomials for Ci := 



and Kac polynomials for 


Ci := 1. We emphasize the difference between the polynomial studied in this 
paper, i.e. the right-hand side of Eq. (10), with the elliptic polynomial: in our 

^ ^ ' are binomial coefficients, not their square root as in the 


case Ci = 


elliptic polynomial. In the former case, v{xy Cv{y) = J2i=i 


d-l (d—1 


xY = 


(1 -I- xy)‘^ and as a result the density function and the expected number of 
real zeros have closed formula, see [15] 


fspit) = , Eep = y/d- 1. 

7r(l -I- E) 

Our case is more challenging, because of the square in the coefficients, 

/d _ ^ 

v{xY'Cv{y) = X)i=i ( j ) ^*2/* is no longer a generating function. Nev¬ 
ertheless, Theorem 2 shows that E(2, d) still has interesting asymptotic be¬ 
haviour as in (18) and (19). 


3.2 Connections to Legendre polynomials and other qualitative properties 

In this section, we first establish a connection between the expected density 
function fd and the well-known Legendre polynomials. Then using the connec¬ 
tion and known properties of Legendre polynomials, we prove some qualitative 
properties of fd and the expected number of equilibria. The main results of 
this section can be summarised as follows. 

(i) Theorem 3 and Theorem 4 derive an expression for the expected density 
fd in terms of the Legendre polynomials. 

(ii) Theorem 5 shows that is an increasing function of d for any given 
t > 0. 

(iii) Corollary 4 proves that and are decreasing functions of d. 

Technically, keys to these theorems are Lemma 3 that connects the Legendre 
polynomials Pd to Md+i in (II), and Lemma 4 showing a reverse Turan’s 
inequality. These lemmas are of interest in their own right. 
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3.2.1 Legendre polynomials 


We recall some relevant details on Legendre polynomials. Legendre polynomi¬ 
als, denoted by Pd{x), are solutions to Legendre’s differential equation 


d 

dx 




-|- d{d + l)Pd{x) — 0, 


(27) 


with initial data Po{x) = 1, Pi{x) = x. The following are some important 
properties of the Legendre polynomials that will be used in the sequel. 

(1) Explicit representation 


= ^E (•)'(*- 

(2) Recursive relation 

{d+ l)Pd+i{x) = {2d+ l)xPd{x) - dPd-i{x). 

(3) First derivative relation 

= xPd{x) - Pd-l{x). 

Example 2 The first few Legendre polynomials are 

Po(a::) = 1, Piix) = x, P2ix) = ^{3x‘^ - 1), Ps^x) = ^{5x^ 

Pa{x) = i(35a;"^ — 30x^ -I- 3). 

8 


(28) 


(29) 


(30) 


3a;), 


The Legendre polynomials were first introduced in 1782 by A. M. Legen¬ 
dre as the coefficients in the expansion of the Newtonian potential [33]. This 
expansion gives the gravitational potential associated to a point mass or the 
Coulomb potential associated to a point charge. In the course of time, Legendre 
polynomials have been widely used in Physics and Engineering. For instance, 
they occur when one solves Laplace’s equation (and related partial differential 
equations) in spherical coordinates [5] . 

Our approach reveals an appealing, and previously unexplored relationship 
between Legendre polynomials and evolutionary game theory. We explore this 
relationship and its applications in the next sections. 


3.2.2 From Legendre polynomials to evolutionary games 

The starting point of our analysis is the following relation between the poly¬ 
nomial Md+i{t) in (11) and the Legendre polynomial Pd. 

Lemma 3 It holds that 

1 + t'^ 


Md+i{t) = {1-ffPd 


1-P 


(31) 
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Proof See Appendix 6.1. 


Note that traditionally in the Legendre polynomials, the arguments are in 
[—1,1]. In this paper, however, the arguments are not in this interval since 
I 1 > 1- Legendre polynomials with arguments greater than unity have 
been used in the literature, for instance in [13, Chapter 2]. We now establish a 
connection between the density fd{t) and the Legendre polynomials. According 
to the second property in Proposition 1, we have 


fd+i{t) 


1 


1 

t V Md+i(t)) 


(32) 


where ' denotes the derivative with respect to t. Using this formula and Lemma 
3, we obtain the following expression of fd+i{t) in terms of Pdit) and its 
derivative. 


Theorem 3 (Expression of the density in terms of the Legrendre 
polynomial and its derivative) The following formula holds 


{27:fd+iit)f 


4d^ 

(l-t2)2 


i6f^ 


(33) 


Proof See Appendix 6.2. 

Corollary 3 As a direct consequence of (33), we obtain the following bound 

for f2,d{t)- 

m < (34) 

In comparison with the estimate (14) obtained in Theorem 1, this inequality 
is weaker for t > 0 since it is of order 0(d— 1). However, it does not blow up 
as t approaches 0. 


We provide another expression of fd+iit) in terms of two consecutive Legendre 
polynomials Pd-i and Pd- In comparison with (33), this formula avoids the 
computations of the derivative of the Legendre polynomial Pd- 

Theorem 4 (Expression of the density fnnction in terms of t-wo Leg¬ 
endre polynomials) It holds that 


(27r/d+i(t))^ 


4d2 d^ll+t^ Pd-l fl+t^'\ 

(1 -t2)2 “ fT [l_i2 “ J 


(35) 


Proof See Appendix 6.3 
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3.2.3 Monotonicity of the densities 


Theorem 4 is crucial for the subsequent qualitative study of the density f 2 ,d(t) 
for varying d. 

Lemma 4 The following inequality holds for all |a;| > 1 

Pd{xf < Pd+i{x)Pd-i{x). (36) 


Proof See Appendix 6.4. 

Note that this inequality is the reverse of the Turan inequality [47] where 
the author considered the case x G [—1,1]- 

Theorem 5 For any given t > 0, is a decreasing function of d. 

Proof We need to prove that 

fd+ljt) ^ fd+2{t) 
d ~ d+1' 

From (35), we have 


47r^ 


1 + t2 
l-t2 


Pd-i 

Pd 


l + t"^ 
1-P 


(37) 


fd+iit) y _ 4 1 

d ) (l-f2)2 

Assume that cc > 1. Then Pd{x) > 0 for all d and from (29), we have 

{d+ l)Pd+i{x) > {2d + l)Pd{x) - dPd-i{x), 

which implies that 

{d+l){Pd+i{x)-Pd{x)) > d{Pd{x)-Pd-i{x)) > ■ ■ ■ > l{Pi{x)-Po{x)) = x-l > 

Therefore Pd+i{x) > Pd{x) for all a; > 1. We first consider the case 0 < t < 1. 
From Lemma 4 for a: = > 1, we have 


p2 


l + t2 


< Pd-1 


1+^2 


P, 


d+1 


1+^2 


or equivalently 

Pd+i\i-ty- Pd \i-ty 


< 1 < 


1 +t2 
1^' 


It follows that 


T + ^2 

1 

P+t2X1 

2 

< 

T + ^2 

Pd 


_l-t2 

Pd 

[l-tyl 

_l-p 

pd+1 



From this inequality and (37) as well as a similar formula for /d+ 2 , we obtain 

fd+ljt) ^ fd+2{t) 
d - d+1 ’ 


which is the claimed property for the case 0 < f < 1. The case t > 1 is then 
followed due to the relation fd{t) = (2 /d (J) ■ 
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Since E{2,d) = fd{t) dt and SE{2,d) = the following statement is 

an obvious consequence of Theorem 5. 

Corollary 4 and are decreasing functions of d. 

The following proposition is a necessary and sufficient condition, which is 
stated in terms of the Legendre polynomials, for f 4 ^i{t) being increasing as a 
function of d. 

Proposition 2 fd+i{t) is an increasing function of d if and only if, for x = 

i+t^ 

l-t2 


{d + lf[Pi+i{x) - P^{x)] ■ [Pi+i{x) + Pj{x) - 2xPd+i{x)Pd{x)] 

+ {2d + l){x^ - l)Pj{x)Pj+,{x) > 0. (38) 

Furthermore, if 


{2d + l)Pi > Pj_, [{2d - + 2P2] (39) 


then (38) holds. 

Proof See Appendix 6.5. 

We numerically verify the inequality (39) in Figure 2; however, it is unclear 
to us how to prove it rigorously. We also recall that it is shown in Theorem 
1 that fd{t) behaves like y/d — 1, which is an increasing function of d, as d 
is sufficiently large. Motivated by these observations, we make the following 
prediction. 

Conjecture 1 For any given t > 0, fd{t) is an increasing function of d. 

We provide further numerical simulation to support this Conjecture by 
directly plotting fd{t) in Figure Ic. 

Remark 2 We recall that in the case n = 2, the variable t is defined hy t = 
where y is the fraction of strategy 1 and 1—y is that of strategy 2. Equivalently, 
y can be expressed in terms of < as y = . Hence one can also transform the 

statements of the theorem above (and later) in terms of y. As has been shown in 
the beginning of Section 2, the transformation from y to t has the advantage 
that it transforms a complex equation (Eq. (8)) to a univariate polynomial 
equation (Eq. (10)). This enables us to exploit many available results and 
techniques from the literature of random polynomial theory. Moreover, from 
the relationship between f{t) and g{y), it follows that all the monotonicity 
properties with respect to d are reserved for 32 ,d (see a numerical illustration 
in Eigure 6). 
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4 General Cases 

In this section, first we prove an estimate for the density /n, 2 (t) similarly as 
in Theorem 1 for a two-player multi-strategy game. The expected number of 
internal equilibria in this case has been computed explicitly in [14, Theorem 
3]. We then conjecture an asymptotic formula for the general case. Finally, we 
provide numerical simulations to support our conjectures, as well as the main 
results in the previous sections. 


4.1 Two-player multi-strategy games 


Theorem 6 (two-player multi-strategy games) 

Assume that {A'' = (/3))j=i^... = 1,... ,n — 1} are independent random 

vectors, then it holds that 


/n, 2 (t) = 7r i/ ti • •-tn-i = 0 and 

/n, 2 (t) < Tr~^r min 1 1 , „ ^ - 1 if h ■ --tn-i 7 ^ 0 . 

V2/ ( n2 ti ■■■ tn-i j 


As a consequence, for any t such that ti • ■ ■ t„_i 7 ^ 0 

lim = 0. 

n->oo TT 2 T 

Proof According to [14, Theorem 3] 

/n,2(t) = TT-^r Q) 


n—1 




k^l 


(40) 

(41) 


(42) 


(43) 


The first assertion is then followed. Now suppose that ti ■ ■ ■ t„_i 7 ^ 0. By the 
Cauchy-Schwartz inequality, which states X]"=i ''' an)" for n all n 

positive numbers ai, • • • , a^, we have 


n-l \ t 

1 + X! tn-l- 

k=l / 


Therefore 


/n,2(t) < TT = A Q) 


n2 ti 


tn.— 1 


/ „-l \ 2 

On the other hand, since ( 1 + X) ) > Ij it follows that 

\ fc=i / 


/n,2(t) < TT 


Therefore, we obtain (41). 
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Note that the expected number of internal equilibria for a two-player multi¬ 
strategy game is given by [14, Theorem 3] 

E{n,2) = ^. 


Table 1 Expected number of internal equilibria for n = 3 and different d, as generated 
from E{3,d) and from averaging over 10® random samples of the system of equation in 2. 
In all cases, results from random samplings are slightly smaller E{3,d). 

_d_2_3_4_ 5 

E{3, d) from Theory 025 057 092 L29 

E{3, d) from Random Sampling 0.249496 0.569169 0.910236 1.28898 


Remark 3 In Theorem 6 , we need an assumption that the random vectors 
{A^ = (/3*)j=i_... ,n-i; * = 1, ■ ■ ■, n—l} are independent. Recalling from Section 
2 that where is the payoff of the 

focal player, with i {1 < i < n) being the strategy of the focal player, and 
ik (with 1 < ik ^ n and 1 < fc < d — 1 ) is the strategy of the player in 
position k. The assumption clearly holds for n = 2. For n > 2, the assumption 
holds only under quite restrictive conditions such as is deterministic 

or a], f. are essentially identical. Note that the assumption is necessary 
to apply [15, Theorem 7.1]. Hence, to remove this assumption, one would 
need to generalize [15, Theorem 7.1]. This is difficult and is still an open 
problem [30,31]. Nevertheless, since the system (3) has not been studied in 
the mathematical literature, it is interesting on its own to investigate the 
number of real zeros of this system under the assumption of independence of 
{H*}. As such, the investigation not only provides new insights into the theory 
of zeros of systems of random polynomials but also gives important hint on 
the complexity of the game theoretical question, i.e. the number of expected 
number of equilibria. In Figure 7 and Table 1, we numerically compare the 
density function 33 ^ and the expected number of equilibria computed 
from the theory with the above assumption and from samplings without the 
assumption. We observe that have the same shape (behaviour) in both 
cases. In addition, the values of E^ d computed from samplings are slightly 
smaller than those computed from the theory. 


We now make a conjecture for the expected density and expected number 
of equilibria in a general multi-player multi-strategy game. 

Conjecture .2 In a multi-player multi-strategy game, it holds that 

hiE{n, d) n—l 


fn,d {d-1)^ 


and 


lim 

d-fcxD ln((i — 1 ) 


(44) 
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We envisage that even a stronger statement holds 

E{n,d) (45) 


This conjecture is motivated from Theorem 2 and a similar result for the sys¬ 
tem of multivariate elliptic polynomial as follows (see also Remark 1). Consider 
a system of n — 1 random polynomials, each of the form 


0<fci <cZ—1, 

ki—d—i 



(46) 


where the coefficients are independent standard multivariate normal 

distribution. Then according to [15, Section 7.2], the expected density and the 
expected number of real zeros are given by 

/^p(t)=^-5r(|)|^^J^, EEP = {d-l)^. (47) 

These formula are direct consequences of Lemma 1 and the fact that v{x)Cv{y) 
is a generating function, which generalises the univariate case 


u(x)^Cu(y) = 


E 

ki—d—l 


d-1 

ki^kj 




2=1 


(l+x-y)‘'-i. 


As mentioned in Remark 1, in our case u(x)^Cz;(y) is no longer a generating 
function due to the square of the multinomial coefficients. Motivated by this 
and Theorem 2 for the univariate case, we expect that the conjecture holds 
for the general case (i.e. d-player n-strategy normal evolutionary games). 


4.2 Numerical results 

In this section, we provide numerical simulations for all the (main) results 
obtained in the previous sections. The following figures are plotted. In the 
following list, (1) to (4) are to illustrate and confirm the analytical results 
obtained in the previous sections. The others, from (5) to (10), are numerical 
simulations for the games with large d and n. They also provide numerical 
support for Conjecture 2. 

(1) fdd) and fd{t) / {d— 1) as functions of t for different values of d, see Figures 
la-b. Figure la illustrates the fourth property in Proposition 1 that fd{t) 
is increasing as a function of t. Figure lb explains the scaling by showing 
that fd{t)/{d — 1) is a bounded and decreasing function. 
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(2) /d(i) and /d(t)/ {d— 1) as functions of d for different values of t, see Figures 
Ic-d. These figures show that fd{t) increases with d while fd{t)/{d— 1) 
decreases, which are in agreement with Conjecture 1 and Theorem 5. 

(3) fd{x) and fd{x)/{d — 1) as functions of the frequency x. 

(4) ^ function of d, see Figure 3a. This figure demonstrates Theo¬ 
rem 2 on the convergence of to 1/2 as d tends to infinity. 

(5) f'i,d{t\,t 2 ) as a function of t = for different values of d, see Figures 

5a-c. 

(6) f 3 ,d{ti,t 2 ) as a function of d for different values of t = see Figure 

4a. 

(1’) / 4 ,d(fi) fa) as a function of t = {ti,t 2 ,t^) for different values of d, see 
Figures 5d-f. 

(8) fi,d{ti,t 2 ,tz) as a function of d for different values of t = (fi,f 2 ,f 3 ), see 
Figure 4b. 

(9) as a function of d, see Figure 3b. 

(10) as a function of d, see Figure 3c. 

In Figures 5, we provide numerical results of /n,d(t) for n = 3 and n = 4. 
We observe that the density function decreases with ti (namely, ti and t 2 for 
n = 3, and ti, t 2 and ts for n = 4) and increases with d. We conjecture that for 
the general d-player n-strategy normal evolutionary game, the density function 
decreases with ti and increases with d. 

Figures 3b and 3c support Conjecture 2. From these two figures, one also 
can see the complexity of the problem when d increases. We are able to run 
simulations for d up to 10000 for n = 2, up to 400 for n = 3 and only up to 
20 for n = 4. 


5 Discussion and outlook 

How do equilibrium points in a general evolutionary game distribute if the pay¬ 
off matrix entries are randomly drawn, and furthermore, how do they behave 
when the numbers of players and strategies change? To address these impor¬ 
tant questions regarding generic properties of general evolutionary games, we 
have analyzed here the density function, fn,d, and the expected number of (sta¬ 
ble) equilibrium points, E{n,d) (respectively, SE{n,d)), in a normal d-player 
n-strategy evolutionary game. We have shown, analytically and using numer¬ 
ical simulations, that f 2 ,dit) monotonically decreases with t while it increases 
with d. The latter implies that, as the number of players in the game increases, 
it is more likely to see an equilibrium at a given point t. We also proved that 
its scaling with respect to the number of players in a game d, , decreases 

with d. More interestingly, we proved that this density function asymptoti¬ 
cally behaves in the same manner as Vd — 1 at any given t (i.e. regardless of 
the equilibrium point). Similar monotonicity behaviors of the density function 
were observed numerically for games with larger numbers of strategies n. Ad¬ 
ditionally, we proved an upper bound for two-player game with an arbitrary 
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Fig. 1 (a) Plot of / 2 ,(i(t) and (b) of 1) as functions of t, for different 

values of d. We observe that both functions decrease with t, which are in accordance with 
Proposition 1. The second function is bounded from above, having maximum at t = 0, which 
agrees with Theorem f. (c) Plot of /2 £j(f) and (d) of /2 d(t)/(d— f) as functions of d, 
for different values of t. We observe that /2 d{t) increases with d, while f 2 ,dit)/i<i ~ 1) 
decreases, which are in agreement with Conjecture f and Theorem 5. All the results were 
obtained numerically using Mathematica. 


number of strategies, n. Briefly, our analysis of the expected density of equi¬ 
librium points has offered clear and fresh understanding about equilibrium 
distribution. Related to this distribution analysis, there have been some works 
analyzing patterns of ESSs and their attainability in concrete games, i.e. with 
a fixed payoff matrix [49,11,8]. Differently, our analysis addresses distribution 
of the general equilibrium points, and for generic games. Note also that those 
works dealt with two-player games while we address here the general case (i.e. 
with arbitrary d). 

Regarding the expected numbers of internal (stable) equilibrium points, 
first of all, as a result of the described monotonicity properties of the density 
function, we established analytically that E{2,d) and SE{2,d) increase with d 
while their scaled forms, and decrease with d. Next, we proved 

a new upper bounds for n = 2, with arbitrary d: E{2, d) ^ y/d — 1 ln((i — 1). 
This upper bound is sharper than the one described in [14] (which is also the 
only known one, to the best of our knowledge). As a consequence, a sharper 
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2 4j^6 8 10 2 ^ X ^ ^ 10 2 4j^68 10 



Fig. 2 Plot of L(i(x) := (2c/+ 1)P^—|^(2f/ — l)Pj_|_j + 2Pjj for different values of d. 
We observe that Ld{x) is always non-negative, thereby supporting Conjecture 1, since the 
inequality L^ix) > 0 is the sufficient condition for the conjecture, as shown in Proposition 
2 



Fig. 3 Plot of as functions of c/, (a) for n = 2, (b) for n = 3 and (c) for 

n = 4. The results confirm the asymptotic behaviour for n = 2 as described in Theorem 
2, and clearly support Conjecture 2 for the general case. All the results were obtained 
numerically using Mathematica. 


upper bound for the number of expected number of stable equilibria can be 
established: SE{2,d) < \\/d — 1 ln(d — 1). More importantly, that allowed us 

to derive close-form limiting behaviors for such numbers: lim 1 

d^oo 

and lim = |. As such, apart from the mathematical elegance of 

d —>-00 ^ ' 


























































24 


Manh Hong Duong, The Anh Han 




d 


d 


Fig. 4 Plot of (a) /3,£i(ilii2) fof different values of (fi,t2) and (b) /4,(i(flif2) for 
different values of (ti,42,43), both as a function of d. We observe that the functions 
increase with d while decrease with ti, with 4 = 1, 2 in panel (a) and 4 = 1, 2, 3 in panel (b). 
All the results were obtained numerically using Mathematica. 


our results, the analysis here has significantly improved existing results on the 
studies of random evolutionary game theory [14,22,16]. Moreover, generalizing 
these formulas for the general case, i.e. when the number of strategies n is also 
arbitrary, our conjecture that, lim is nicely corroborated by 

numerical simulations. Related to this analysis, there have been some works 
analyzing ESS in random evolutionary games [19,29,25]. In both [19] and [29], 
the authors focused on the asymptotic behaviour of ESS with large support 
sizes, i.e. considering also equilibrium points which are not internal, in random 
pairwise games. In a similar context, in [25], the authors studied ESS but 
with support size of two, showing the asymptotic behaviors of such ESS when 
the number of strategies varies. Differently from all these works (which dealt 
exclusively with two-player interactions), our analysis copes with multi-player 
games. Hence, our results have led to further understanding with respect to the 
asymptotic behaviour of the expected number of stable equilibria for multi¬ 
player random games. 

Last but not least, as the density functions we analyzed here are closely re¬ 
lated to Legendre polynomials, and actually, they are of the same form as the 
Bernstein polynomials, we have made a clear contribution to the longstanding 
theory of random polynomials [6, 15,45,38] (see again discussion in introduc¬ 
tion) . We have derived asymptotic behaviors for the expected real zeros Eg of a 
random Bernstein polynomial, which, to our knowledge, had not been provided 
before. Note that the asymptotic behaviors and close forms of the expected 
real zeros of some other well-known polynomials have been derived by other 
authors. For instance, the Weyl polynomials for Ui '■= the elliptic polynomi¬ 
als or binomial polynomials with := 



; and Kac polynomials with 
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Fig. 5 Plot of (a-c) as a function of (£i,t 2 ); and (d-f) / 4 ,d(ti,i 2 ) as a 

function of (ti,t 2 ), for different values of ts. We observe that both functions increase 
with d while decrease with ti, with i = 1, 2 in panels (a)-(c) and i = 1, 2, 3 in panels (d)-(f). 
In panels (a) and (d), d = 2; in panels (b) and (e), d = 5; in panels (c) and (f), d = 10. In 
panels (d)-(f), the surfaces, from bottom to top, correspond to £3 = 0.1, 0.3, 0.5, 0.7 and 
10, respectively. All the results were obtained numerically using Mathematica. 


Gi := 1. The difference between the polynomial studied here (i.e. Bernstein) 
with the elliptic case: ai = ^ are binomial coefficients, not their square 

root. In the elliptic case, v{x)^Cv{y) = ^ ^*2/* = (1 + 

and as a result, E{2, d) = y/d— 1; While in our case, because of the square in 

/d _ ^ 

the coefficients, v{xYCv{y) = ^ j x''y'’, is no longer a generating 

function. Indeed, due to this difficulty, the analysis of the Bernstein polyno¬ 
mials is still rather limited, see for instance [3] for a detailed discussion. 

In short, our analysis has provided new understanding about the generic 
behaviors of equilibrium points in a general evolutionary game, namely, how 
they distribute and change in number when the number of players and that of 
the strategies in the game, are magnified. 


6 Appendix 

Detailed proofs of some lemmas and theorems in the previous sections are 
presented in this appendix. 
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(a) (b) 



Fig. 6 (a) Plot of g 2 ,d{y) and (b) of g 2 d{y)/{<^~ 1) as functions of y (note that 
t = for different values of d. We observe that both functions are symmetric about 

the line y = 1/2, in accordance with Lemma 2. The first function increases with d while the 
second one decreases. All the results were obtained numerically using Mathematica. 






Fig. 7 (Top row) Plot of the density function g^ ^(^ 15 ^ 2 ) for different values of d (here 
yi and y 2 are fractions of strategy 1 and 2, respectively, and 1 — yi — y 2 is that of the third 
strategy). (Bottom row) Probability histograms plots of results from solving the system of 
equations 2 with n = 3 and for different d, where the payoff entries a’s are (independently) 
sampled from Gaussian distribution with mean 0 and standard deviation of 0.5 (10® samples 
in total). The probability histograms have similar shapes to g 3 ,d{yi-,y 2 ) (see also Table 1). 
All the results were obtained numerically using Mathematica. 


6.1 Proof of Lemma 3 


This relation has appeared in [18, Exercise 101, Chapter 5]. For the readers’ 
convenience, we provide a proof here. 
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Proof (Proof of Lemma 3) Let x = from (28) we have 



By taking q = , we obtain (31). 


6.2 Proof of Theorem 3 


Proof (Proof of Theorem 3 ) 

By taking the derivative of both sides in (31), we obtain 




1 + 


+ Tt{i - p'^ 


i + t^ 

l-t2 


-2 1 d +4t(i - t^y-^p'^ 


1+^2 


It follows that 

M'+i(t) -2td At P'/l + t2\ 

Md+yt) + 


Now we compute the expression inside the square-root of the right-hand side 
of (32). We have 


-2yd Ay p(/i + y\ 

^ {i-yyp,yi-y)^ 

and 

f.K+iy)\_ ^td 8yi + y)pwi + y\ ley p(fp,-{p(^y /i + y\ 

\Md+yt)) (l_i2)2+ (l_i2)3 (^I_i2y + (l_i2)4 p2 - y J ' 
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Substituting this expression into (33), we get 


(27r/d+i(<))^ 


Ad 8{1+ t^) P' + 16t^ P'lPd-{P',? 

Pd {i-py Pj 


(48) 

According to (27), the Legendre polynomial Pd satisfies the following equation 
for all a; € R 


l + t^\ 
1-p) ' 


-2xPd{x) + (1 - x^)Pd{x) = -d{d + l)Pd(a;)- 
As a consequence, we obtain 




1 


2 a; 5 M-d(d+l) ) ■ 


Pd{x) 1 - a;2 Pd{x) 

Substituting this expression into (48) with x = , we get 

4d 


{2T:fd+i(t)Y = - 


(l-t2)2 

16t^ 


%{l + e)Pd + 

(l-<2)3p^ 


1 


a-'^)Mi-(si)^ 


+ fi + e 


Ad^ 


16P 


p> 

(1_P)2 {l_t2Y\Yd 

which is the claimed relation (33). 


1 - t2 P, V 1 - p 


l+P 


l-P 


-d{d+l) 



6.3 Proof of Theorem 4 


Proof (Proof of Theorem 4) Using the following relation of the Legendre poly¬ 
nomials for all a; S R 


= ^ 2 ! 1 i^Pdix) - Pd-i{x )), 

we get 

P'dix) ^ d / _ Pd-i{x) \ 
Pd{x) a;2 - 1 \ Pd{x) J ' 

In particular, taking x = we obtain 


( 49 ) 


Ek(\±f\ 

Pd \l-P) 


d 


'l + t^ 



- 1 


l-P 


Pd-i A + 

Pd 


dii-t^Y 


+ 

l-t2 


Pd-1 

Pd [i-ty 


AP 
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Substituting this expression into (33), we achieve 


(27r/d+i(<))^ = 


AdP (f\l+e Pd-i fl + i 

(l-t2)2 “ ^ [l-t2 “ J 


which is (35). 


6.4 Proof of Lemma 4 


Proof (Proof of Lemma 4) This lemma follows directly from [12] (Theorem 
2.1) where the authors proved that 


Pdix)^ - Pd+iix)Pd-iix) = 


1 — 


d-1 


did+l) 1 ^ z + 1 

\ 1=1 1=1 J=1 


which is negative for all |a;| > 1. 


6.5 Proof of Proposition 2 


Proof (Proof of Proposition 2) We will prove that 

/J+2(i)-/i+i(0>0 ^ (38). 

From (35), we have 

4^^(/i+2(t)-/J+i(t)) = 


(50) 


4(d+l)2 

(d+l)2 

'l + P 

Pd 

(l-t2)2 

P 

_l-P 

Pd+i 


4^2 d2 


(1 -t2)2 

4(2d+l) 1 

(l-t2)2 “ p 


1+^2 

l-t2 


Pd-1 

Pd 


1 +t2 

l-t2 


1+^2 

l-t2 


id+l)Pi-dPd-lPd+l 


PdPd+1 


1 + ^2 




d + dPd-iPd+i f I +1 


PdPd+l 

Therefore fdit) is increasing as a function of d if and only if 


2 

l-P 


1+^2 (d + 1)P^ — dPd-iPd+i f 1 + t 

1 — ^2 PdPd+1 

4(2d+l)t2 


2 \ 1 


l-t2 


{2d + 1) 


1 + t2 (d + ^)Pd T dPd-iPd+i /1 + t 


l-t2 


PdPd+1 


^-2\l 


1-P 


< 


{1-py 


We re-write the expression above using the variable x, using the relation x'^ — 
1 = as follows 


(l-t2)2> 

(d -|- l)Pj — dPd-iPd+i 
PdPd+l 


{x) 


{2d + l)x - id+l)Pi + dPd-iPd+i 

PdPd+l 


(a;) 

(51) 


< (2d-hl)(a:2-l). 
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We now simplify this expression using the recursion relation of the Legendre 
polynomials, i.e. dPd-i = (2d + l)xPd — (d + l)Pd+i- Namely, we have 


xPdPd+i — (d + l)Pd + dPd-iPd+i — xPdPd+i — (d + l)Pd + [(2d + l)xPd — (d + l)Pd+i]Pd+i 

= (d + 1 ) [—+ 2xPdPd+i] 


and 


{2d + l)xPdPd+i — (d + l)Pj — dPd-iPd+i 

= {2d + l)xPdPd+i - (d + 1)P2 _ [(2d + l)xPd - (d + l)Pd+i]Pd+i 
= (d+l)(Pi+i-Pj). 

Substituting these calculations into (51) we obtain (38). 

To prove the second assertion of Proposition 2, we proceed as follows. Let 

Hd+i = {d+ if [Pd+i{x) + Pj{x) - 2xPd+i{x)Pd{x)\ 

= [(2d + l)xPd{x) - dPd-i{x)f + (d + lfPj{x) 

- 2x{d + 1) [(2d + l)xPd{x) - dPd-i{x)] Pd{x) 

= [{2d + Ifx"^ + (d + If - 2{d + l)(2d + l)xf Pj{x) 

— [2d(2d + l)a; — 2d(d + l)a;] Pd{x)Pd-i{x) + d^P'^_fx) 

= [d^ + (2d + 1)(1 - x^)] Pd{x) + d‘^Pd_i{x) - 2xd^Pd{x)Pd-i{x) 
= [Pd{x) + Pj-i{x) - 2xPd{x)Pd-i{x)] + (2d + 1)(1 - x^)Pj{x) 
= Hd + {2d+l){l-x^)Pj{x). 


Hence, the expression in (38) can be simplified as follows 


Hd+i[Pl+fx) - Pj{x)] + (2d + l){x^ - l)Pj{x)Pj^fx) 

= [Hd + {2d + 1)(1 - x^)Pj{x)][Pj+fx) - Pj{x)] + (2d + l){x^ - l)Pj{x)Pj+fx) 
= Hd[Pj^fx) - Pj{x)] + (2d + l)(x2 - l)Pi(x) 

Pj+fx)-pj{x) 


Pii^) - Pi-ii^) 


[Hd {Pj{x) - Pj_fx)) + (2d - l){x^ - l)Pj{x)Pj_fx) ■ Q] , 


where Q 


(2d+l)Pl{Pl-Pl_^) 

{2d-l)Pl_,(Pl^,-Pl)- 


Suppose that (39) is true, i.e.. 


(2d + fPt > [(2d - l)Pj+i + 2Pj] . 
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This implies that Q > 1 for all x and d. Then it follows that 
Hd+i[P!+,{x) - P^{x)] + {2d + l){x^ - l)Pi{x)P!^,{x) 

> [Hd {p!{x) - PlAx)) + {2d - l){x^ - l)P^{x)Pl,{x)] 


^Pl^{x)-P!{x) P]{x)-Pl,{x) 

- Pj{x) - pj_,{x) Pj_,{x) - Pl,{x) 


X {Pl,{x) - Pj_^{x)) + {2d-3){x^ - l)P^_,{x)Pl^{x)] 


> 


^ n X [^1 - Poi ^))+• 

(52) 

By definition of Hd^ we have 

Hi = Pi{x) + Pq{x) — 2xPi{x)Po{x) = x^ + \ — 2x^ = 1 — x^ . 
Substituting this into (52), we obtain 

Hd+i[Pj+i{x) - Pj{x)] + {2d + l){x^ - l)Pj{x)Pj+i{x) 

■ Pli{x)-Pf{x) 


> (a^"-i)n 




Pf{x)-Pli{x) 


= Pl^i{x) - P2 {x) > 0, 
i.e., the condition (38) is satisfied. 


Acknowledgements We would like to thank the anonymous referee for his/her useful 
suggestions that helped us to improve the presentation of the manuscript. 


References 

1. N. H. Abel. Memoire sur les equations algebriques, ou I’on demontre I’impossibilite de 
la resolution de I’equation generale du cinquieme degre. Abel’s Ouvres, (1):28—33, 1824. 

2. L. Altenberg. Proof of the Feldman-Karlin conjecture on the maximum number of 
equilibria in an evolutionary system. Theoretical Population Biology^ 77(4):263 - 269, 
2010. 

3. D. Armentano and J. Dedieu. A note about the average number of real roots of a 
bernstein polynomial system. Journal of Complexity, 25(4):339-342, 2009. 

4. R. Axelrod. The Evolution of Cooperation. Basic Books, New York, 1984. 

5. S.S. Bayin. Mathematical Methods in Science and Engineering. Wiley, 2006. 

6. A. Bloch and G. Polya. On the Roots of Certain Algebraic Equations. Proc. London 
Math. Soc., S2-33(l):102, 1932. 

7. M. Broom, C. Cannings, and G. T. Vickers. On the number of local maxima of a 
constrained quadratic form. Proc. R. Soc. Lond. A, 443:573—584, 1993. 

8. M. Broom, C. Cannings, and G. T. Vickers. Sequential methods for generating patterns 
of ess’s. Journal of Mathematical Biology, 32(6):597—615, 1994. 












32 


Manh Hong Duong, The Anh Han 


9. M. Broom, C. Cannings, and G.T. Vickers. Multi-player matrix games. Bull. Math. 
Biol, 59(5):931-952, 1997. 

10. M. Broom and J. Rychtar. Game-theoretical models in biology. CRC Press, 2013. 

11. C. Cannings and G. T. Vickers. Patterns of ess’s ii. Journal of theoretical biology, 
132(4) :409-420, 1988. 

12. E. Constantinescu. On the inequality of p. turan for legendre polynomials. J. Inequal, 
in Pure and Appl. Math, 6(2), 2005. 

13. L. J. Curtis. Atomic Structure and Lifetimes. Cambridge University Press, 2003. 
Cambridge Books Online. 

14. M. H. Duong and T. A. Han. On the expected number of equilibria in a multi-player 
multi-strategy evolutionary game. Dynamic Games and Applications, pages 1-23, 2015. 

15. A. Edelman and E. Kostlan. How many zeros of a random polynomial are real? Bull. 
Amer. Math. Soc. (N.S.), 32(l);l-37, 1995. 

16. C. S. Gokhale and A. Traulsen. Evolutionary games in the multiverse. Proc. Natl. Acad. 
Sci. U.S.A., 107(12):5500-5504, 2010. 

17. C. S. Gokhale and A. Traulsen. Evolutionary multiplayer games. Dynamic Games and 
Applications, 4(4);468—488, 2014. 

18. R. L. Graham, D. E. Knuth, and O. Patashnik. Goncrete Mathematics: A Foundation 
for Computer Science. Addison-Wesley Longman Publishing Co., Inc., Boston, MA, 
USA, 2nd edition, 1994. 

19. J. Haigh. The distribution of evolutionarily stable strategies. Journal of applied prob¬ 
ability, pages 233—246, 1988. 

20. J. Haigh. How large is the support of an ess? Journal of Applied Probability, pages 
164-170, 1989. 

21. T. A. Han, L. M. Pereira, and T. Lenaerts. Avoiding or Restricting Defectors in Public 
Goods Games? Journal of the Royal Society Interface, 12(103):20141203, 2015. 

22. T. A. Han, A. Traulsen, and C. S. Gokhale. On equilibrium properties of evolution¬ 
ary multi-player games with random payoff matrices. Theoretical Population Biology, 
81(4):264 - 272, 2012. 

23. The Anh Han. Intention Recognition, Commitments and Their Roles in the Evolution 
of Cooperation: From Artificial Intelligence Techniques to Evolutionary Game Theory 
Models, volume 9. Springer SAPERE series, 2013. 

24. G. Hardin. The tragedy of the commons. Science, 162:1243-1248, 1968. 

25. S. Hart, Y. Rinott, and B. Weiss. Evolutionarily stable strategies of random games, and 
the vertices of random polygons. The Annals of Applied Probability, pages 259-287, 
2008. 

26. C. Hauert, S. De Monte, J. Hofbauer, and K. Sigmund. Replicator dynamics for optional 
public good games. J. Theor. Biol, 218:187-194, 2002. 

27. J. Hofbauer and K. Sigmund. Evolutionary Games and Population Dynamics. Cam¬ 
bridge University Press, Cambridge, 1998. 

28. S. Karlin. The number of stable equilibria for the classical one-locus multiallele slection 
model. J. Math. Biology, 9:189-192, 1980. 

29. S. C. Kontogiannis and P. G. Spirakis. On the support size of stable strategies in random 
games. Theoretical Computer Science, 410(8):933—942, 2009. 

30. E. Kostlan. On the distribution of roots of random polynomials. In From Topology 
to Computation: Proceedings of the Smalefest (Berkeley, CA, 1990), pages 419-431. 
Springer, New York, 1993. 

31. E. Kostlan. On the expected number of real roots of a system of random polynomial 
equations. In Foundations of computational mathematics (Hong Kong, 2000), pages 
149-188. World Sci. Publ., River Edge, NJ, 2002. 

32. E. Kowalski. Bernstein polynomials and brownian motion. Amer. Math. Monthly, 
113(10):865-886, 2006. 

33. M. Legendre. Recherches sur I’attraction des sphrodes homognes. Mmoires de Math- 
matiques et de Physique, prsents VAcadmie Royale des Sciences, par divers savans, et 
lus dans ses Assembles, Tome X, 1785. 

34. J. Maynard-Smith. Evolution and the Theory of Games. Cambridge University Press, 
Cambridge, 1982. 

35. J. Maynard Smith and G. R. Price. The logic of animal conflict. Nature, 246:15—18, 
1973. 



Title Suppressed Due to Excessive Length 


33 


36. A. McLennan. The expected number of nash equilibria of a normal form game. Econo- 
metrica, 73(1):141-174, 2005. 

37. A. McLennan and J. Berg. Asymptotic expected number of nash equilibria of two-player 
normal form games. Games and Economic Behavior^ 51(2):264-295, 2005. 

38. H. Nguyen, O. Nguyen, and V. Vu. On the number of real roots of random polynomials. 
Communications in Contemporary Mathematics, to appear, 2015. 

39. M. A. Nowak. Evolutionary Dynamics. Harvard University Press, Cambridge, MA, 
2006. 

40. J. M. Pacheco, F. C. Santos, M. O. Souza, and B. Skyrms. Evolutionary dynamics of 
collective action in n-person stag hunt dilemmas. Proc. R. Soc. B, 276:315-321, 2009. 

41. M. Perc and A. Szolnoki. Coevolutionary games — a mini review. Biosystems, 99(2): 109 
- 125, 2010. 

42. S. Petrone. Random bernstein polynomials. Scandinavian Journal of Statistics, 
26(3):373-393, 1999. 

43. F. C. Santos, M. D. Santos, and J. M. Pacheco. Social diversity promotes the emergence 
of cooperation in public goods games. Nature, 454:213-216, 2008. 

44. P. Schuster and K. Sigmund. Replicator dynamics. J. Theo. Biol., 100:533-538, 1983. 

45. T. Tao and V. Vu. Local universality of zeroes of random polynomials. Int Math Res 
Notices, 2014. 

46. P. D. Taylor and L. Jonker. Evolutionary stable strategies and game dynamics. Math. 
BioscL, 40:145-156, 1978. 

47. P. Turan. On the zeros of the polynomials of Legendre. Casopis Pest. Mat. Fys., 
75:113-122, 1950. 

48. G. T. Vickers and C. Cannings. On the number of equilibria in a one-locus, multi allelic 
system. J. theor. Biol., 131:273-277, 1988. 

49. G. T. Vickers and C. Cannings. Patterns of ess’s i. Journal of theoretical biology, 
132(4) :387-408, 1988. 

50. E. C. Zeeman. Population dynamics from game theory, in Lecture Notes in Mathemat¬ 
ics, 819:471-497, 1980. 



